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ABSTRACT 


Equations are developed which give the pressure profile, the forces 
and torques on a disk pendulum by means of point source wave theory from 
acoustics. The pressure, force and torque equations for an unbaffled 
disk are developed. These equations are then used to calculate the 
apparent mass and apparent inertia for the pendulum. 
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other sources on a disk 

value of r for source at edge of disk 
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upper limits of integration on r 
impedance 

normal to boundary 

angle of swing of pendulum 

wavelength 

ambient air density 
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Subscripts : 

a apparent 

eg center of gravity 

i incident at a differential area 

s area of pressure source 

s' area where pressure is applied 

Dots over symbols denote differentiation with respect 
to time. 

Hats over symbols denote maximum values. Roman numeral, 
suboripts I and II denote pressures due to transla- 
tion and rotation, respectively. Arabic subscripts 
denote successive terms in the total pressure. 



INTRODUCTION 


An expression for the total pressure on the surface of a rigid 
circular disk vibrating in a baffle vas worked out by Rayleigh 
(reference 1, Vol. II page 162). However, Rayleigh did not give an 
expression for the pressure distribution across the disk. This was 
done subsequently by McLachlen (reference 2, p. 1012) by means of an 
expansion using hypergeometric functions. Such a pressure distribution 
is necessary if correct moments are to be calculated for a disk in 
rotation about a diameter as frequently occurs in experimental work. 

The work of these authors can be extended readily, in the 
special case of low frequencies, to include the pressure on a disk 
pendulum. The pendulum can be considered as an unbaffled disk in 
rotation and translation. The net affect of the absence of a baffle 
is to reduce the pressure on both faces since air is free to travel 
around the disk to the opposite face as the pendulum swings. 

The pressure acting on the face of a disk represents _n inertial 
force if dissapative effects are neglected. Hence the pendulum in 
motion appears more massive than when not in motion. The result is 
that the period of the pendulum will be different at different air 
densities. Physically the pendulum drags a certain amount of air with 
it as it swings. This air is referred to as attached mass, entrained 
mass, apparent mass, or virtual mass by various authors. For a 
pendulum, rotational effects must also be considered; so apparent 
inertias are also present as shown in sketch (a). 
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Pendulum 

motion 


Translation 4* 


Rotation + Pendulum 
length 


Sketch (a) 


It will he the objective of this paper to show that the theoreti- 
cal acoustic pressure profile can be investigated by breaking the 
equations of motion into two special cases: (a) a pure translation 

and (b) a pure rotation about a center line of the disk. A combination 
of these effects, along with consideration of the pendulum support 
distance, affords a description of the radiation and attached mass 
by means of the principle by superposition and hence represents the 
total effect. 

A comparison will be made of these results with the theoretical 
results obtained from classical hydrodynamics and from experiment. 
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Theoretical Mod»l 

In the present study two models can be considered as technically 
feasible. One of these consists of ignoring the absence of a baffle 
and thus treating the pendulum as a di ributed monopole sound source 
at low frequency. In this case the Green's function becomes the free 
space Green's function for outgoing waves with infinite boundary conditions 
in the positive half space. Under these conditions the pressure is finite at 
the edge of the disk. The second model consists of assuming that the 
pendulum is a distributed acoustic dipole. This leads to a somewhat 
more complicated model than a monopole disk since under these condi- 
tions a term must be included in the Green's function for pressure which 
serves tc. force the free spece term to zero at the boundary. By 
symmetry it also follows that the pressure must be zero at all points 
in the plane of the disk beyond the boundary. 

In the present 1 .per a dipole model will be used to calculate the 
attached mass of the disk. 

The present paper is similar in some respects to a number of 
papers already in the literature. However, several of significant 
differences are discernable. M. StraBSberg (reference 3, p. 520) 
considered the radiation field only at large distances from the 
oscillating body. The result is that the force may be computed 
correctly but no detailed description of the pressure distribution is 
possible, although an approximate value can be obtained by setting 
F * ma. Bouwkomp (reference b ) and Wiener (reference 5) treated the 
sound fields diffracted around various obstacles. As such their 



papers are valuable in that they complement the present work. 

However, in both papers, fixed integration limits result in a more 
complicated sumnation formulation which oust satisfy matching condi- 
tions at the boundary and, as such, these papers are more like pertur- 
bation problems than is the present paper. Crane (reference 6) takes 
an approach which is closest to that of the present work in that a 
free space potential field is considered along with an auxiliary 
field to satisfy the boundary conditions. However, crane assumes 
zero potential from the edge of the disk to infinity whereas in the 
present paper this condition arises out of symmetry considerations 
and the pressure is guaranteed zero only at the boundary. This results 
in Crane following Bouwkomp's paper by insisting on a fixed coordinate center 
for integration. The resu.lt is that in references 4, 5 and 6 it becomes 
necessary to describe the pressure from a coordinate system which is 
centered at the center of the disk; and hence these authors are forced 
to use sunmations of Bessel and Legendre functions. In the present 
paper the pressure is calculated in a coordinate system located in 
each instance at ehe assumed receptor point. While circular symmetry 
is sacrificed by this approach the algebra becomes such thr.t the only 
expansions required arc those of the ordinary exponential. It should 
be noted that in none of the above references is rotational moti;. 
about a disk diameter considered, hovever, Mangulis (reference 7) 


does consider this for a baffled disk. 



MATHEMATICAL MODEL 


Monopole Disk Theory 

Consider, first, a rigid disk Titrating in an infinite baffle 
vith simple harmonic notion. The total force is determined by inte- 
grating the pressure over the surface of the disk. The pressure at 
sene infinitisimal area ds' produced by the oscillation of some 
nearby area ds is given by equation (l) 


d Pc 


3 p 0 ck a «J <ut - ,tr) 

kin- 


ds 


( 1 ) 


vhere r is the distance between ds and ds', u is the Martini 
speed of the oscillator, k the wave n<mber, c the speed of sound, 
Pq the ambient air density, and S the surface area of the disk. 
Thus the total pressure acting on one face can be obtained by 
integrating over the surface of the disk treating each infinitis i mal 
area ds' as a point source (reference 8, p. 178). 

The total reaction force acting on the disk is the integral of 
he pressure over the surface. 



Substituting equation (l) into equation (2) yields 
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Figure 1.- Coordinates employed in describing the acoustic pressure at 
ds due to a source at ds' for a disk of radius a. 
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F = 



, ~ J(wt - kr) 
cku e" 

5irr 


ds 


ds' 


(3) 


The coordinate system for performing the indicated integration over 
the front s urface of the disk is ahovn in figure 1 (references 8 , 9 » 
and 10 ). 

From figure 1, it can be seen that 


2 

r -a cos 4 + [a 
max — 


sin 


1/2 


(U) 


It can also be seen by inspection of special cases, th«*t the + sign 
is the correct selection for all quadrants given the convention for 
4 as shown in figure 1. Hence for all situations 


max 


* o cos 


4 + [a 2 - o 2 sic 2 4 ] 


1/2 


(5) 


With this expression in the limit on r the pressure at ds 
due to an sources on the front surface of the disk becomes, in polar 
coordinates centered at ds 


ncos 

u u 


4 + [a 2 - a 2 sin 2 4] 


1/2 


. /> - kr) 

J p Q cku e 

5tt 


dr d$ 


( 6 ) 


This equation can be written in a form where the limit is slightly 


modified as 
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-n 


a[— ccs ♦ ♦ (1 - 


a 2 2 1/2 

<!> ♦» jp 0 c«. 


J(«t - kr) 

1 - dr d$ 

( 7 ) 


The Green's function for one side of & disk radiating to the half 
space in an infinite baffle (reference 11) is 

J(wt - kr) 

G , *-S — (8) 

monopole 2^r 

Dipole Disk Theory 

If the effect of the back face of the disk is also considered the 
pressure at the edge, i.e. when r * 2a cos must be zero since the 
wave produced by one face is exactly 180° out of phase with the other 
and tb geometries are identical. Thus we specify a general Green's 
function by 

J((»it - kr) j(wt - k[ba cos $ - r]) 
dipole ” Virr TO" cos $ - rj ™ 

( 9 ) 

It can be seen that this choice of a general Green's function does in 
fact produce the desired pressure at the edge. 

There are several important points to be considered in relation 
to this Green's function. 

1. The function G, . , satisfies the characteristics of a 

dipole 

Green's function; i.e., the first term on the right contains the 
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essential singularity vh-le the first and second term together satisfy the 


boundary conditions and do not contain any singularities over the region. 

2. With this choice of a Green's function the pressure can be 

determined by the application of Green's theorem since on the boundary both 

G and the outward normal gradient - 5 — are zero vhile the pressure is 
0 dn Q 

continuous. The line integrals over the boundary therefore vanish; and 

only the surface integral contributes to the pressure. 

3. The physical significance of the Green's function is the 

following: The first term on the right of equation (9) is a monopole 

source located at any point on the front face of the disk. The 
second tern is a monopole located along r on either the front or 
back face of the disk having strength Just sufficient to cancel the 
source on the front face at the opposite boundary. Some examples are 
shown in figure 2 . 

The pressure on the front face of the disk is now given by 
integrating G over the surface of the front face 


P = 


i P 0 cku 
Ttt 


2x a 

n 


„ ,2 , l/ 2 r 

a[— cos $ ♦ (l - (— ) sin‘ $) ]i 

A & ft 


ri kr 


-Jk[Ua cos $ - 


[Ua cos $ 


- r] 

s — r d r 

- r] 


d* 


( 11 ) 


By a change of variable let 


x 



( 12 ) 


also let 
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( 13 ) 


for convenience. These changes reduce the foregoing equation to the 
simpler form 


J P 0 cu 

P = <*> 


2 2 

2n m cos $ + (l - m sin 


n 



e -j(ka)[l» cos ♦ - x] 
" [4 cos <J> - x] 


x dx d<t> 


(14) 


Since the second term in the integrand must be expanded in a 
Taylor ' s series before integration, there is no advantage to be 
gained in not expanding both terms. The expansions are, for (ka) « 1, 
i.e. very long wavelength 


e -j(ka)x 

x 


and 
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Source on front face, 
image on back face 
and receptor on 
boundary 


Source and image on opposite 
faces of boundary, receptor 
on opposite boundary 


Image on 
front face 



Incomplete 

cancellation 



Image on 
back face 

Incomplete 

cancellation 


Source and image on Source on front face, image 

front face and receptor on back face and receptor 

half way between on front face 


Figure 2.- Source and image arrangements for a dipole disk 
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Hence, to sufficient accuracy 


, „ A Jut 
J P 0 cu e J 

p ~ Uir 


2 2 

2ir a cos $ + (l - m sin $) 


1/2 


(ks> J J | jv COB l - Si * ii i l ~ (l * 00* ♦ - SxTj Xix 4* 


J p Q cu e' 

CtT 


Jut 


0 0 

2tt m cos $ ♦ (l - m^ sin^ $)^^ 

(ka) j J J ■ [bx cos $ - 8 cos^ $]x dx d$ 


0 0 


The second integrand on the right above is real and corresponds to 
energy being radiated away from the disk. Both terms in the first 
integrand are imaginary and represent pressure due to attached mass. 

It can be seen that for ka « 1 the second attached mass term is 
negligibly small in comparison to the first. Hovever, as this term is 
larger in magnitude than the first radiation term it will be carried 
in the present analysis for the sake of completeness. The integrals 
will be taken one at a time. 


First Attached Mass Term 
The first term in equation (17) can be written 


P 1 * 


.1 p Q cu e 


0 ut 


m cos $ + (l - m^ sin^ 


(ka) 


2jr m coj 

JJfr 

•'o •'o 


ax a» 

cos $ - x) 


(19) 


The iirst integral of this can be found directly from tables and 


is given by 
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p l * 


J P Q cu e 


Jut 


Un 


(ka) 




cos 4> + 2(1 - m 2 sin 2 $)*^ 2 


+ cos $ An 


1 _ a cos 6 ♦ (l - m 2 sin 2 
U cos <(> 


d(j» (19) 


Since an integral of the natural logarithm term in equation (19) is 
difficult it becomes expedient to expand this term to make the 
expression more tractable. An appropriate and highly accurate 
expansion is the standard form 


An(l - x) = - [x + — + 




( 20 ) 


[x 2 < 1 and x - -l] 


This expansion is accurate over most of the surface of the disk 
as can be seen from examination of figure 3. Also since x is a 
function of cos $ the integral of this expansion exists and is 
finite. Then 


J p cu e'^ 

W 


(ka) 


r- 


cos + 2(1 - m 2 sin 2 4») 1 ^ 2 


, /. 2.2, »l/2 [m cos d> + (1 - m 2 sin 2 d*) 1 ^ 2 ] 

- a cos <t> - (1 - m sin $) - J 2(i"eo« » ) 


_ fm cos d> + (l - m 2 sin 2 ift) 1 ^ 2 ] 
2(4 cos $)^ 


+ higher order terms} d<J> (21) 



Raising terms to the indicated powers and carrying out the integration 
gives 


p l = 


J P Q CU e 


Jut 


hrr 


(ha) 


E(f , . <^>K<f. .,] 


( 22 ) 


Where K(^, m) is the complete elliptic integral of the first kind 
and where E(^, m) is the complete elliptic integral of the second 
kind. Now making use of the standard approximate expansions of the 
respective elliptic integrals 


2 4 

K(^, m) 5 f (1 + f” + + * * * ^ 


(23) 


and 


E(~, m) S £ (1 - J- - + . . • ) 


(24) 


and carrying terms to the fourth power in m 


J P Q cu e 


Jut 


(ka)[l - ■^■ni-^m^ + ^g-m^ - m* + higher order terms] 


(25) 


A plot of this pressure profile is given in figure 4 across a radius 

of the disk. In addition, a pressure profile for the hack surface of 

o 

the disk is identical to that of the front surface except for a 180 




16 


phase shift. The total pressure difference is thus obtained by 
do x the pressure on one side. 


Second Attached Mass Term 
The second term in equation (17) can be written 


J p Q cu e^ wt 
P 2 (to) 


2tr a cos $ + (l - m 2 sin 2 <p)*^ 


IS 

0 0 


(ka) 2 


[4 cos $ - 2x] x dx d. 


(26) 


The first integral of this is 


. „ ^ jurt 

J p_ cu ~ 

-a ’ —Hi < to > 3 


fr 

I l cob <j> [m 


cos <|> + (l - m 2 sin 2 <fc)^ 2 ]^ 


[m cos <fr + (l - m 2 sin 2 &)^ 2 1 3 *^ 


J 


d<t> 


(27) 


Raising terms to the indicated powers and carrying out the indicated 
integrations over <(> gives 


J p Q cfi e Jwt . 2 k 2 it 

P 2 * — *— (ka ) 3 {[fj - fj(l + m 2 ) B(§, m) 


| (1 - n 2 ) E(| s m) - - |](1 - m 2 ) K(|, m) (28; 


It can be seen that this term is smaller than the preceeding term 

o 

by a factor on the order of (ka) where ka « 1 at low frequencies. 
Hence this part of the pressure term is completely negligible in 
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ir relation to the first attached boss tern. If the elliptic integral 
expansions, equations (23) and (2U), are employed as before and the 
results carried up to the fourth power in m the pressure becomes 


J P n c “ 




Po “ 


(ka)' 




m ~k 


m - 


1 3 
lS“ 


1 

65: 


* J (29) 


A scaled plot of equation (29) is shown in figure 5. 


First Non-Cancelling Radiation Term 
Strictly speaking the firs* radiation terms from the front and 
back faces of the disk cancel each other. In the present analysis 
this cancellation was done in the expansion of equation (lU) . Hence 
the normal monpole radiation, proportional to (ka), is not present 
for a dipole. The first non-cancelling radiation term is therefore 


P 3 = " 


p Q cu e 


Jut 


3* 


(ka) 


2ir m cod $ + (1 - m P sin^ <J>)^ P 

'/J 1 


g 

[x cos $ - 2 cos $] x dx d$ 


(30) 


0 0 


The integral of this expression over x is 


a Juxt 

°° I, l 1 * 

P 3 j; (ka) 


2if r , 

C J (m cos d * (1 - m 1 

H 


2 .in 2 4) 1 ' 2 ! 

■ — . - * ■ cos <|> 


- [m cos $ + (1 - v? oin^ <(>)^^] cos^ $ ) u<J> 


(31) 


Squaring and cubing as indicated, and int' - - iting again over 2tt gives 



Elliptic integral 
soln. is singular 
at x ■ 0 
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aanssdJd 


Figure 5.- Acoustic pressure due to the second attached mass 
r id ns . 


20 


p cu e Jwt u 1 

P 3 * 3 (k*)** [1 - a «■ | 


(32) 


Elliptic integrals do not appear in the real terns. 

It can he seen that for ka « 1 chis tern is even smaller in 
magnitude than the previous two terns. However, as this is the first 
term which radiates energy away from the disk, it is important from 
a theoretical standpoint. A scaled plot of (32) is given in figure 6 . 


Total Pressure for a Disk in Translation 


The total pressure for a disk in translation without a baffle is 

a 


given by p ♦ p ♦ p and is explicitly, with — for m 

X c j & 

(l - -* 


a Jwt/iiiai (i 9 . T.. 12 . ii- -"LiiL \ E (i £ 

k* \ 12 / M 2 * a 

^ L 


p = p Q cu e' x 


) 


- V 


(1 - (J) 


12 


J w£ £) + f(2a . ]l 

/ K *2* a' w U 3o 9 


)d ♦ <f> > «f. 


- 1 (i - (f) 2 > e(|. f) - (§ - !><i - tf) 2 ) «f. -)) 


♦ U - (2) ♦ i (|) 2 ]j 


( 33 ) 


ka « 1 


This same equation with the ellipx ' c integrals expanded and 


carried to the fourth power in (~) is 

& 
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P * % * u - K> - $ <f ) 2 * fe (|) 3 * £ <fn 

* J ^ [- i ♦ <*> - £ (f) 2 - ^ (2) 3 ♦ & (f) 1 *] 


+ 




(3**/ 


Force on a Disk in Translation. 

The force on one face of the disk is obtained by integrating the 
pressure over the surface of the disk. Using equation (3U) 



With the indicated operations carried out, the force is 
F ® p Q cu e Jwt a 2 {j[l.068l (ka) + 0.2962 (ka) 3 ] + [0.6109 (ka) 11 ]} 

(36) 

Since the two faces of the disk are out of phase the total force 
is twice this, or 

F = p n cu e* 50 * a 2 {j [2.1362 (ka) + 0.592 1 * (ka) 3 ] + [1.2218 (ka) U ]} 

° (37) 
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Acoustic Impedance 
The radiation impedance is given by 


Z 

r 


F 


i Jwt 
u e" 


( 38 ) 


Hence the radiation impedance of the disk iB given by 


Z r = p Q ca 2 {j[2.1362 (ka) + 0.5924 (ka) 3 ] + [1.2218 (ka) 1 *]} (39) 


as a point of comparison, the radiation impedance of a simple dipole 
is veil knovn. In terms of the present notation this impedance is 
given (ref. 11, p. 317) by 


Z « p ca 2 {j[2.09U3 (ka) + 1.0472 (ka) 3 ) + [1.047 (ka) 1 *]} (40) 

r u 


It can be seen that there is reasonable agreement between the disk, 
vh>ch is a distributed dipole, and a simple dipole at low frequency 
as would be expected, with the diBk slightly the higher of the two. 

Apparent Mass 

The apparent mass is given by the non-radiating, or imaginary 
part of the reaction to the force. 


F 



v) 


(41) 



Or, if the apparent mass is assumed constant 


F = 


n 

a dt 


In the case of a disk in translation 


(U2) 


dv 

dt 


- J wue 


Jut 


(1*3) 


Therefore, using equations (1*2), (1*3), and the first imaginary term 
of (37) 


m = 2.1362 p n a 3 (1*1*, 

a 0 

Pressure on a Disk Pendulum 

That the effective mass and inertia of a disk pendulum are greater 
in air than in avacuum can be readily demonstrated experimentally. 
However, in the past, these effects have generally been measured 
mechanically (ref. 12 and 13) and theoretical treatment has been by 
means of classical hydrodynamics. The result led to some difficulty in 
visualizing the relative part played by rotational terms and transla- 
tional terms, especially when the pendulum is of intermediate length 
so that rotation and translation contribute significantly to the total 
motion. An alternative approach is taken in the present paper; viz., 
to consider the pendulum as an acoustic source, albeit of extremely 


low frequency. 
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Consider a pendulum which consists of a disk suspended by & 
massless support and swinging in still air about an axis in the plane 
of the disk, so that the air impinges on the flat face of the disk. 

If the motion of the pendulum is undamped, as it would be if the 
pendulum were suspended in a vacuum, then the equation of motion for 
small amplitude oscillation is 


16 + m g l 8 = 0 
s c.g. 


( 45 ) 


The period of oscillation is given by the well known formula 


P = 2tt 


v= 


8 8 ^C.g. 


( 46 ) 


with circular angular velocity 


AS* 


( 47 ) 


In contrast, the damping moment acting on a pendulum oscillating 

in still air is proportional to the drag (fig. 7) which, in turn, is 

proportional to the square of the velocity. In thi3 case, the drag 

1 2 

force is given ’ y the product of the aerodynamic pressure x pV , 
the surface area of the pendulum on which this pressure acts 3, and 
the drag coefficient of the pendulum C^. The damping moment is the 
product of the drag force and the moment am. Hence, for small 
amplitudes, the equation of motion with drag added is given by 
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“ * | “ v l v l s C D l aero * “= « ‘eg. 6 * 0 < k6 ) 

The moment of inertia I includes the structural and apparent moments 
cf inertia. However, the restoring moment is related to the 
structural mass m g only. The apparent mass produces no restoring 
force, since the buoyant force balances the gravity force. 

The velocity term in equation (1*8) is given by 


V 


aero 


Therefore, equation (48) can be written as 


(49) 


« * | s °E «l®l ♦ “ s * l c .g. 8 * 0 <50> 


It should be noted that this equation differs from the more commonly 
encountered equation of damped oscillatory motion in which the damping 
is proportional to the first power of the velocity. 

Equation (50) is a nonlinear, second-order, differential 
equation which doe3 not heve e known solution. However, the behavior 
of this equation has been examined numerically for a wide range of 
constant damping coefficients. A typical solution is shown in 
figure 8 which was taken from reference 13. It can be seen that the 
period is nearly constant, is independent of the damping, and is equal 
to the period which would be predicted on the assumption that no 





Figure 8.- Typical time histories (eq. ,(h)) for damped and undamped oscillations. Init xl 
conditions: * 0.2 rad, = 0 rad/sec. 
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damping is present at all. The same result was obtained in the 

approximate solution of equation (30) given in appendix B. These 

results have also been verified experimentally by the author (not 

published). The results, therefore, indicate that, if the damping is 

proportional to the square of the velocity, the period is essentially 

the same as that for the undamped situation. Hence, the circular 

velocity is given to sufficient accuracy >y equation (1*7). The value 

of this is that the period in air gives a measure of the attached mass. 

In this equation m , g, and 1 are constants. For purposes of 

s c .g. 

this analysis, it will be assumed that I depends only on the mass and 
geometric parameters cf the pendulum and on the density of the attached 
air. Hence I can also be considered constant and w is constant. The 
disk of the pendulum mny, therefore, be considered as a distributed 
acoustic source oscillating at angular rate w. Two possible types of 
motion will be present, viz. translation and rotation. Both types of 
motion will produce pressure effects on the face of the disk. These 
pressure effects can be broken up into forces and torques in the 
following manner: 

Pressure Across the Disk of the Pendulum 
The maximum speed u of any infinitesimal point on the surface 
of the disk is determined by its distance from the pendulum support 

A 

point. Thus if 0 is the maximum angular speed of the pendulum, 
i.e.; as it passes through bottom dead center 
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U = u - b )0 (51) 

c.g. 

where b is the distance of the point where the pressure is applied 
from a horizontal line through the center of gravity of the disk, 
(figure 9). 

From figure 9 

b * r sin (<p - ip) + o sin ip 

Then (51) becomes 

A /V A 

• • ♦ 

u - l 0 - a sin ip 0 - r sin (4> - ij>)0 (51a) 

C • g* 

This equation can be inserted in equation (11) so that 



0 0 


This can be broken into two integrals, 


-Jk[U a cos <}> - r] 


n " ' — T ~ — T~ r dr d$ 
[4 a cos $ - r] 
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Center line 
of pendulum 


Figure 9*- Coordinates used to specify speed as a function of distance 
from the center line of pendulum. 
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1/2 


SL “ ' q ' — ~ — t * 

[“ f i p Q ck[£ c g< e-o sin ^ Sle^* 

-Jkr 

e 

e “Jk[k a cos <P - r' 

J J ^ 

r 

[U a cos <p - r] 

_ 


r dr 


0 0 

2 1/2 

2 ir a[^ cos <f> + (l - f— ) sin^ $) ] 

| ^ p Q ck[r sint^jeje^ 

0 0 


IF 


e -jfcr e -Jk[4 a cos $-r] 


[U a cos $ - r] 


r dr d4> 
(53) 


The first of these two integrals is seen to he the sane as equation (11) 


except for the substitution of the term [ & 6 - a sin i|> 0] in place 

C *^t 

of the constant u. Hence the integral of this expression is merely 


equation (33) or its expanded form (3*0 with [£ 0-0 sin \|> 0] 

C *g * 

replacing u. That is, part of the pendulum pressure profile has the 
same form as that of the simple disk developed previously. If (53) 
is written in the convenient form 


p i = P I + p II 


(5U) 


with Pj analogous to (33) 


p i * p o ol •'•o.g.® - 0 sin * u - r <f> - k ‘f* 

* fe<f> 3 ♦ 1 ♦ ^ t- 1 ♦ <f> - t<r )3 - is<f> 3 * 



( 55 ) 
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and 


P II 


2rr a[— cos $ + (1 - (— ) sin 2 $) ] 

a a 


f f J 

c— ' 

e -Jkr 

e ~Jk[U a cos $-r] 

3 3 55 

r 

[h a cos $-r] 


0 0 


sin($-^ ) dr d$ 
(56) 


Again there is no particular advantage to be gained in integrating 
the first term in the square bracket directly so both terms are 
expanded in a Taylor’s series. First change the variables by letting 

x * - (12) 


and 


0 

m * — 
a 


(13) 


as was done previously. Then equation (56) becomes 


v* ' 3Mt , ■ 


2 2 

2 ii m cos <|» + (l - m sin 


1/2 


II 


u sin(<HI') 
0 0 


g -J(ka)x e -j(ka)[U cos 4>-x) 


Pt cos <(> - xj 


x 2 dx d$ 


(5T) 


The term in brackets can now be expanded as before with the result 


that 
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2 2 

2i r m cos $ + (l - m sin 4>) 


1/2 


P II = 


J P Q c0 e 
l*ir 


Jut 


ka 


| | sin(4H>) 

0 0 






+ 1 (4 cos 0 - 2x) 


x dx d<fr 


2 2 

2n m cos + + (1 - m sin <J>) 


i/2 


Jut 


Uir 


J P 0 cB e— " f f ( } 3 2 2 

ka I pi^^-^) *•* • [1* x cos $ - 9 cos $]x dx d$ 


0 0 


(58) 


Again there are two imaginary terms in lower powers of (ka) before 
reaching the first real radiation term. Again, the integrals will he 
taken one at a time. Designating these integrals as p^, and pg 
respectively, one has the following. 


First Attached Mass Rotation Term 
The first term is 


2 2 

2it m cos $ + (l - m sin $) 


1/2 


J p ° c§ ^ 
h m R to 


fj 


sin(^) 


(4 cos $ - 2x) 
x(4 cos $ - x) 


x d dx d$ (59) 


0 0 


The first integral of this gives 
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J p c6 e J f 2 x 

^ ka sin(^) x m +kx n cos cos *) In 1 - 


where is the upper limit 


2 2 J 

x = m cos i ♦ (1 • i sin $) 
m 


as was done in equation (20) £n(l - x) is expanded 


2 3 k 5 

0/1 \ r.X.X.X.X. 

ln(l - x) = - LX + — + — 4 ~ + T - + 


[x^ < 1 and x = -l] 


This expansion is accurate over most of the surface of the disk 
as can be seen from examination of figure 3. Hence to three term 


accuracy 


J p ee e J r 2 

ka I sin(4> - ty) [x^ + kx m cos <j> 


- {k cos 4>)‘ 


f x m 4 

W cos ♦> * 2(1. co, 


2 3 

cos $) cos $) 


+ higher order terms \ ] d$ 
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A cancellation in the x tern yields 

z& 



If x is written explicitly in terns of $ and use is made of the 
m 

identity 


sin($ — ) — sin $ cos 4* - cos sin i> 


(65) 


the indicated integration over 2ir can be performed. The result is 



If the elliptic integrals are expanded by (23) and (2k) and terms up 
to the fourth power in m are retained the above equation becomes 


^ m 2 - m^ + m^Jsin 4» (67) 


Second Attached Mass Rotation Term 


The second term is 
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P 5 = 


J P q c 8 e' 


Jwt 


2 2 

2ir m cos $ + (1 - m sin $) 


1/2 


Uir 


ka 


J jsin(M») 

0 0 




(U cos $-2x)| z dx d<t> 

(68) 


The first integral of this expression is 


j p.c0 e JUt 


2n 


P 5 = - 


8ir 


(ka) 3 a J sin(<HMj cos Q x 3 - j x^jd$ 


(69) 


where 


x 

m 


= m cos 


♦ + (i - 


2 . 2 

a sin 


♦ > 


1/2 


as before. 

Then writing (69) explicitly in terms of $ rad making use of 
(65), (69) can be integrated. The result is 

a j 

p = - — ° - ^ - e [(12 + 6m - 112m 2 + 8fcm 3 - £8m 4 + 6m 5 )E(£, m) 

5 ^ U5 m 2 2 


(12 + 6m - 88m 2 + 36m 3 + 76 a* - h2a ?) K(^, m)] sin \j> (70) 


as before, the elliptic integrals can be expanded and terms up to m 
retained to yield 


* i p n c0 e Jut (ka) 3 a [7 


s 


( 71 ) 
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First Radiation Term in Rotation (That Does Net Cancel) 
The first non-cancelling radiation tent is 


p 6 = “ 


p ° c ° eJut (k i 1 * 
r (ka) a 

3^ 


2 2 

2tt m cos $ + (1 - m sin <J>) 


1/2 


if 


o 2 2 

sin($-4>)[x J cos <p - 2x cos <|>]dx d<}> 


0 0 

The first integral of this is 

2ir 


(72) 


p 6 


P c c0 e Jwt 
3v 


r 

(kaj^a I sin($ - cos <P - cos^ $]d<|> (73) 


with 


2 2 

x ffi * m cos $ + (l - m sin $) 


1/2 


As before, making use of (65) and writing (73) explicitly in terms of 
<J>, the equation can be integrated. The result is 


p 6 


• Ju>t 

(k.) 1 ** $ . | „ . ♦ I m^lsin * 


(7U) 


as usual, the radiation term does not contain elliptic integrals. 

Pressure Distribution of a Pendulum 

The final result of the foregoing analysis is the total pressure 

distribution at low frequency of a disk pendulum of arbitrary size 

and length. Cne has merely to add up the pressures being careful to 

c 

substitute [A - a sin in place of u in p- , p 5 , and p 0 . 
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The pressure is 


J p c[Jt - o sin \J >]0 
_ 0 c.g. 

p _ ^ 


(1 - (f) 2 ) _ _ 

12 K( 2’ a' 


(ka) 

Jut 




J P n c8 e* 

+ 0»)a 


-{ 


2 3 

12 - fc(£) . 12(2.) + U(J) 

a a a 

9 (f) 


( 


12 - 7 (r) + i 2 (r> - (r 


9 (f) 


£ ) \ 
— ) 


I 

•) «f- f> 


sin <|> 


j P n c[fc - a sin i |/]0 . _ 

* g e^tka) 3 


(-^ -£>(1* (f) Mb f) 


- 5 * 1 - <f> )«?• f> ' 


L 3( £> 


2 l 9 


U - (f) )K(f , f) 


Jut 

'0"\~ (ha) J a 


j P-C0 e' , “ w /._n3_ _ - 2 „ 3 _ U 


2n 


_ [(12 + 6(£) - 112(7) + 8M7) - 28(£> 

U5(-)2 a a a a 

a 


+ 6(f) 5 )E(J, 7) - (12 + 6(7) - 88 (7) 2 + 36 ( 2 -)" + l6{ Z) k 
a a a a a a 

- W( f, 5 ,K(f,f, 1 sin * . .XWix - (f)< t>f) £ i 


P„o 0 e Jut 


3 (*»> k 4 - f(f) ♦ (f) * |(f) 1 sin 


( 75 ) 


If this same equation is written in a fom vith the elliptic integrals 
expanded the result is 
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P * 


3 PgCfft. gt ~ q s * n 409 


(^)tl - £<f> - ^ ♦ &f) 3 


/w 

5 t a \\ i P ° C ^ eiUt (v i r 1 i ~( a i 1 t a \ 2 1 / a i 3 

' 256^ 1 + 5 {ka)a[ - I + 2( I } “ 


i n k J P n c[£ ~ - ° sin 

+ ^a 3 3 sin ♦ + '^~2 


e Ju,t 2 

(*0 3 [- (f) - £<f) 


- &f> * k<?> i - * “o' 5 «i- I ♦ - 3$‘ 


* fcf> 3 - 3fe«H - * * « • > W U 


PftC8 e 


Jut 


^ --— (ka) U a[£ - §<f) + (?f + ^Hsin * 


-<? + ?I> 3 + 3 


2 , _ 3 

a- 3'a 3 


(76) 


Vertical Pressure Profiles 

Equation (76) contains a special case of particular interest. 

If l is assumed to be zero this equation reduces to a disk 
c.g. 

suspended along its centerline. Under these conditions the motion 
is a pure rotation. Of course, a real disk under these conditions 
would have no restoring moment, hut this could be overcome by properly 
weighting the lower edge of the disk. The special case of (76) which 
results is 

P-C0 e JU)t f r. , „ _ 2 Q _ 3 

^ * 3< to > {h * \ V - h V ‘ 32 ‘a 1 + 2st < X ) J 

* 3(ka)3 {j 3*^.) *{$ -f^f* 

- te) k - i(f) ♦ |(f) 3 j] 


sin ij> 


(77) 
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The pressure has a simple sine dependence on hut a somewhat 
more complicated dependence on radius. The terms of (77) are plotted 
in figures 10, 11, and 12 in scaled form with ip » ^ . Hence these 
plots represent vertical pressure profiles across the pendulum. Since 
no translation is present the pressure given by (27) appears along with 
its mirror image for ip - - ^ as a torque on the disc which gives rise 
to an apparent inertia. As usual, it can he seen that at low 
frequencies, ka « 1, the higher order terms are insignificant in 
relation to the first term. 


Force on the Pendulum 

The effect the pressure has on the pendulum can he broken up 
into pure forces and torques. The force is obtained by applying 
equation (76) for the pressure in equation 2. In this case 
2 tt a 
F 

0 0 


A 

H j p n c[Jl 0 - a sin i|/] , . , „ _ 2 , 3 

fr — ■> - I<f) - §<?> - 


. J V s 1 


' 25? 1 * 3 * 2( ? ' K’ ' ^ ,8in * 


j Pn c I 0 - O sin I|>]e 

+ 2 2 z&i 


Jut 


1 Mf> - 


J Pr.ce e 


jut 


r -'- u *> 3 a %< 2 - |< f ) 2 - 2 (|) 3 - ^ fn . m * 


+ p c[4 0 - a sin (|<3e^ Wt (ka) , *[l - (-) + hfj-) ] 

0 c « g * a 2 a 


„ • m 

P-C0 e 


(ka) U a[£ - §<f) +(ff ♦ |(|)"]sin i|»}o da d* 


3 


(78) 


6i at 
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Figure 10.- Vertical pressure across a disk in rotation due to the first attached inertia 
term as a function of scaled disk radius. 
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psi 90 deg 

Second order term Exact elliptic 

_______ integral pressure 
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Since the order of integration is immaterial in this instance, 
it is advantageous to integrate over first. In this case all 
terms containing sin drop out \ and the force equation becomes 

F * - V K.J f f 1 - E<f> - Wi* * T5<f> 3 - 

0 ^ 

♦ [. i . (2) . . i^i) 3 < Jj<an 


♦ Wl [i . (2.) ♦ A#)) a do 
3 a 2 a f 


(79) 


This is seen to be identical to equation (35) except that u is 

/v 

replaced by £. § . Thus with the second integration carried out, 

C • g • 

the force on one face is found to be 


F = — p_c£ 9 e Jwt a 2 {j[l.068l(ka) + 0.2962(ka) 3 ] [O.Sl&fra) 1 *]) 

0 C.g. 


( 80 ) 


The force on both faces is twice this, or 


F = -p cl I e Jwt a 2 { < j[2.1362(ka) + 0.592l*(ka) 3 ] + [l.22l8(ka) U ]} 
u c .g. 


(81) 


Attached Mass of the Pendulum 


Since the acceleration of the pendulum is given by 



the attached mass is found by dividing the first term of (8l) by (82) 


m * 2.1362 p-a 3 
a 0 


(83) 


This is exactly tue same as was found for a simple dish in translation. 


Torque on the Pendulum Disk 

The torque equation is obtained by integrating the pressure 
equation over the disk vith the appropriate moment arm. If the moment 
arm is selected from the centerline of the disk 

2n a 


t 



p(o sin i|))o da d^ 


( 84 ) 


0 0 


or, applying equation (76) for pressure 



T = - 



j °n c t* " 0 8in *10 «^ Mt 


J ^ _ w BAfti y JV V « ^ 1 <) _ *- T i- 

0 C.R* vr, l/0\ ll/0\ X /O \ 

(ka ti - ^ r ) - 55< ; ) ♦ jj(j) 


. ^2, k * iV ^ 

256V 


5 0»W- ^ ♦ 2(f) - frf) 2 - t<|) 3 ♦ * 


J p ft c[£ - o sin i|»]0 e^ Wt (ka)^ . . 2 , _ 3 ,_ 

* 8 — i 1 ♦ <r> - -istf * 


- 1 Po'9 *!- k * r<f> - §<I> + Va> - ife*!' lEiD * 

♦ P ° C[t C-& , 3 ' ° Si ° 1111 J ,1*^,; . (2, * 

P-.C0 e^ Wt h! r, a 2.3 

+ -S-J (ka) 4 a[i - |(f) + (f) + f<f) ]sin A (a sin *)o do d* 

(85) 


As before, since the order of integration is immaterial , the best 

policy is to integrate over ip first. In this instance terms 

containing £ drop out of the integral and one has 

c »g* 



U8 


a 

- , l I . „ -& .M (fc«> t „3 lff“ 11 o 5 1 0° 5 o 7 , 

{-jp 0 c»# Vto 

J & & & 


J D 0 c9 e 1 


Jut 


J P Q c0 e' 


jut 


/, , r . 2U 3 1 o 1 * lo 5 , 1 cr*\ 

o»)a t- — + — -it— -t~ + gnr 1 

a a a 

\3 r 1 -3 . o” Iff 5 1 o° . 1 ff 7 ! 

(k *> ! -3° t r-5"2-i? - 3*?irT ) 

a a a 


. . i J *A,. .3 r 1 .2 » l o 3 1 0* „ 2_ ff 5 1 a 6 , 

" 3 v 6 e (!») » t- 5 o ♦rr-f~ + 7?i-i8oT> 

a 


p.c0 e 


jut 


3 J» 


X k rl 2 3 o.o ^ 1 a\x 

(ka) a^o - f — ♦ -g + J — 3> ** 

a a 


(86) 


If this equation is integrated again and simplified the torque due to 
the pressure acting on one face of the disk is 


t = Trp Q c0 e jwt ka 5 {j[0.03T2 - 0.10U6(ka) 2 ] + [0.0565(ka) 3 ]} (87) 


The total torque due to the pressure on both faces is twice thiB, or 


t = Jrp Q c0 e^ wt ka 5 {j[O.OTM» - 0.2092(ka) 2 ) + [0.1130(ka) 3 ]} (88) 


Apparent Inertia 


The apparent inertia is given by the reaction to the torque. 

In the case of a simple disk in translation there is no torque and. 
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therefore, nc apparent inertia. However, in the case of a pendulum 


T.I f 
a at 


(09) 


and since 


dt J 


(90) 


The real part of I , to first order, is given by 

a 


I a = P Q t| it a 3 ] a 2 [ 0.0279) 


(91) 


or 


I * 0.1169 p_a 3 
a u 


The apparent inertia due to both faces is twice this, or 


1 * P ft [x * a 3 ] a 2 [0.0558) 
a u 3 


(92) 


or 


I = 0.2338 P n a' 
a u 


(93) 
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Figure 13.- Comparison of the first attached mass pressure term for a dipole with the free 
space Green's function for the same disk as a function of scaled disk radius. 


DISCUSSIOM OP RESULTS 


The results of the present investigation can be compared with 
know results from similar investigations in several instances. First, 
the impedance from a distributed dipole is somewhat larger than that 
which would be calculated from a point dipole of equivalent strength 
as was shown in equations ( 38 ) and (39) • The reason for this is that 
a distributed disk, dipole acts to some extent as its own baffle. That 
is, pressure waves originating in the center of the disk have to travel 
farther to get around the edge than is the case for a point dipole. In 
doing this the sound is attenuated by — and thus cancellation be- 
tween the out of phase components is less pronounced. If this is the 
case, then an upper limiting case should be available in the face Space 
Green's function model since this is equivalent to an infinite baffle 
with no phase cancellation. A comparison of the dominant first order 
pressure along a radius between the free space model given in Appendix 
A and the dipole model is given in figure 13. It can be seen that the 
pressure due to the dipole drops off more rapidly all along the 
radius. If more terms had been carried in the expansions the dipole 
pressure would be exactly zero at the edge of the disk. That is, of 
course, not true for the free space model since zero boundary condi- 
tions were not specified in the Green's function. A similar comparison 
of pressures is given in figure lh for pure rotation. Here the pressure 
builds more rapidly and to a higher value before falling off near the 
edge. The result is that a dipole model predicts less torque since the 
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Figure lU.- Comparison of the vertical pressures between a dipole model and the correspond! 
free space term as a function of scaled disk radius. 
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maximum pressure is applied near the center of the disk with a 
shorter effective lever arm. 

The results of the present analysis can also he compared with 
experiment. Table I (modified from reference 13) gives experimental 
values of (ka) for the pendulum experiment done in this reference. It 
is apparent that for all measurements ka < .006. Hence the pendulum 
of this experiment should be described acoustically to sufficient 
accuracy by the first order terms developed in the present paper. No 
direct acoustic pressures were measured in this experiment. However, 
the attached mass and inertia were calculated. These parameters are 
given in Table II. In addition accepted values of these parameters 
from hydrodynamics are given. In this latter instance (reference lh) 
the attached mass is calculated from the kinetic energy and free 
stream velocity of air flowing past a disk for a non-viscous fluid. 

The conclusion to be drawn from this table is that the dipole 
model gives the lowest attached mass value of all models. 

Since the experimental values are higher than the values of any 
of the theoretical models it seems likely that additional air is 
dragged along with the pendulum which remains unexplained. Some of 
this mass is likely to lie beyond the edge of the disk where only the 
model given by Lamb is rigorously applicable. Seme support for this 
idea is to be found in that the free field model of Appendix A 6u.so 
gives higher values, and it was pointed out previously that this model 
does not presuppose zero p* assure at the edge of the disk. An 
additional physical effect not accounted for by the present theory is 



the possibility of vake formation behind the disk as it strings. 

Since the period of the pendulum is not measurably different for 
amplitudes between approximately thirty degrees and the limits of 
visual perception this is not considered significant. However, the 
possibility of vortex shedding in a vake is not to be completely 
discounted. For the same reason direct friction effects on the 
pendulum are not thought to be significant. However, it is quite 
possible that attenuation of the pressure wave as it travels around 
the edge of the disk occurs. This could be of two forms. Either 
the wave shape could spread so that peak amplitude is lowered, or the 
wave could fail to diffract around the edge of the disk for reasons 
which are not at present explained. It is not thought that node 
formation is a factor in the present situation since the effect would 
be to lower further the already low attached mass values (reference 2). 



CONCLUSIONS 


It is feasible to describe the motion of a disk pendulum by 
breaking the motion up into pure rotation and pure translation and 
treating the surface as a distributed acoustic oscillator operating 
at very lov frequency. However* the dipole model analysis of the pre- 
sent paper gives lover values for attached mass and attached inertia 
at lov frequency then either the hydrodynamic model of Lamb or 
previous experiment undertaken by the author. The conclusion is 
reached that, although this model is reasonable from a qualitative 
standpoint, it is only marginally adequate in its present form to 
explain the experimental results. The model is also not in satisfactory 
numerical agreement with hydrodynamic theory. 


55 



SUGGESTIONS FOR SUBSEQUENT RESEARCH 


Since the present analysis gives a theoretical value for pressure 
at any point on a disk pendulum it would be useful to know actual 
experimental relative pressure profiles across a large disk pendulum. 

The scaled nature of the theory suggests that it would not be necessary 
to obtain absolute numerical results since pressure relative to pressure 
at the center of the disk vould be satisfactory; thus a simple experi- 
ment could be devised. It would also be useful to have wind tunnel 
pressure profile data for a simple disk in a very low speed air stream 
at various angles of attack. Wind tunnel data at the extremely low 
speeds suggested have not been readily available in the past. This 
information should be useful in determining the effects of wake and 
vortex action if present and should also give an idea of the air flow 
pattern across the leading and trailing faces of the disk. 

It is also suggested that a study be made of the effect of the 
test facility walls an the necessity for working in a finite chamber, 
the wavelength is one to three kilometers in reference 13 » may impose 
additional boundary conditions not taken into account in the present 
study. In addition this problem is amendable to scaling as has been 
done in the present paper by means of the factor ka. Thus it should 
be possible to obtain results a higher frequencies for shorter wave- 
lengths as, for example, by loudspeaker experiments. 
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TABLE I.- SUMMARY OF EXPERIMENTAL DATA USING TWO OATA- REJECT ION CRITERIA 
[Disk radius = 0.9525 meter) 


X 

meters 

Speed of 
sound 
c, 

m/sec 

a 

ka 

Mass, m, 
*9 

Density, p, 
kgrm 5 

w m 

20 data-rejection criterion 
<a> 

1.5o data-rejection criterion 
<bl 


Standard deviation, 
o. sec 

Mean period, 

■_ < 

Standard deviation, i 
o. sec 

1098.46 

ESI 



4.42 

0.00161 


3220 

0.000 

3.220 

0.000 

1215.32 


.00572 

.00545 

4.42 

.0339 

weighted slug 


.000 

3.230 

.000 

HE53 

341.395 


.00532 

4.42 

.287 

of mass 

3.297 

.005 

3297 

.005 

1142.48 

339.218 


.00523 

4.42 

.609 

0.850 kg 

3.373 

.010 

3.375 

.010 

1167.89 

339.598 


.00512 

4.42 

.875 


3.436 

.009 

3.440 

MVS 1 


338.429 



4.42 

1.24 

c.g. 

3317 

.018 

3317 

.018 

1521.35 

351.790 

pw! 

0.00393 

4.21 

0.00150 


4.323 

0.005 

4.323 

0.005 

1468.03 

335.327 

UPift? 

moixm 

4.21 

.0349 


4.375 

.005 

4.375 

.005 


341.395 


.00373 

4.21 

287 


4.690 

.000 

4190 

.000 


336.739 


.00351 

4.21 

.618 


5.060 

.000 

$.060 


1795.20 

337.483 

Bmw 

.00333 

4.21 

.886 

437 

5.307 

.010 

5.307 

.010 

1886.84 

338.429 

Ki 

.00317 

4.21 

1.24 

437 

5.582 

016 

5582 

.016 j 

1712.C3 

335.520 



mm 

0.00166 

mam 

5.1® 

0.000 

5.100 

B 

1745.33 

337.755 

.00360 

.00343 

ran 

.0044 


5.168 

.008 

5.165 

.005 

1875.58 

339.624 

.00335 

.00319 

4.21 

.290 


5324 

.008 

5.524 

.005 


341.185 

besi 

.00292 

4.21 

.602 


$.990 

.018 

5,987 


2137.13 

340.572 

.00294 

.00280 

4.21 

.870 


6285 

.015 

6.290 


2284.79 

341.192 

.00275 

.00262 

4.21 

1.22 


6.691 

.010 

6.691 

BBB 

1945.26 


0.00323 

0.00308 

4.21 

0.00161 

8.04 

5.705 

0.005 

BS 

0.005 

1969.65 


.00319 

.00304 

4.21 

.0338 

8.04 

5.775 

.016 


.012 

2115.55 


.00297 

.00253 

4.21 

.289 

8.04 

6208 

,016 

6J2M 

.011 

\fmm 

341.468 

.00274 

.00261 

4.21 

.601 

8.04 

6.705 

.037 

6.705 

.037 

2425.94 

339.596 

.00259 

.00247 

4.21 

.875 

8.04 

7.140 

.016 

7.140 

.016 

2575.08 

338.429 

.00244 

.00232 

4.21 

124 

8.04 

7.594 

.019 

7.600 

.000 

2un.it 

340.689 

0.00303 


mm 

0 00)61 

9.14 

6,088 

0.008 

6.089 

0.007 

2094.40 

340.237 

.00300 


Hfl 

.0339 

9.14 

6.159 

.008 

6.159 

.008 

2268.30 

341.395 

.00277 

.00264 

4.21 

287 

9.14 

6140 

.016 


.013 

2260.14 

339.041 

.00278 

.00265 


291 

9.14 

6168 

.020 

KB 

.005 

2444.82 

IMI 

.00257 

.00245 

ran 

.609 

9.14 

7208 

.018 

BEX 

.010 

2435.34 


.00258 

.00246 

421 

.614 

9.14 

7.195 

.018 

BEX 

.000 

2585.67 

339.792 

.00243 

.00231 

K1 

174 

9.14 

7396 

.044 

BEX 

.010 

2575.08 

338.247 

.00244 

E33 

flt'* $6 

182 

9,14 

7102 

.016 


.016 

270.75 

338.429 

.00229 

.00218 

m 

124 

9.14 

8.099 

.015 

■9 

.015 


fa) Data are rejected it the absolute value of the deviation exceeds two times the standard deviation. 
CdOata are rejected if the absolute value of the deviation exceeds 1.5 times the standard deviation. 



















































TABLE II— SUMMARY OF APPARENT MASS AND APPARENT INERTIA 


RESULTS OF SEVERAL INVESTIGATIONS 



m 

I 


a 

a 

Lamb (Hydrodynamic ) 

2.66T P Q a 3 

0.355 P Q a 5 

Dunning (Dipole) 

2.136 p Q a 3 

0.23k P Q a 5 

Dunning (Free space) 

2.699 P Q a 3 

0.2T9 P Q a 5 

Dunning f Experimental) 

(2.9269 * 0.1195)P 0 a 3 

(0.3933 + 0.0127)p Q a 5 
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APPENDIX A 


Monopole Disk Model 

It is useful for comparison purposes to determine the results for 
a pendulum using only the free space Green's function. Under those 
conditions the pendulum would be so arranged that pressure on one face 
has no effect on the other face. Such would be the case, for instance, 
for an idealized pendulum with the faces isolated from each other by 
an infinite baffle. Since the translational part of the normal pendu- 
lum analysis corresponds so closely with a distributed dipole disk, it 
stands to reason that a baffled pendulum would correspond closely to 
the well known example cf a distributed nonopole as, for instance, a 
loudspeaker in a baffle. The pressure distribution for this configuration 
is well known (reference 2). 

A starting point i* to be had with equation 7* 


2 1/2 
2ir a[— cos $ + (l - (— ) sin $) ] 


- 1 ! ! 


a p 0 ck a 


"tr 


dr d$ 


0 0 


(7) 


Since the first radiation term is not self-cancelling in this case it 
is sufficient to retain only the first order power of k. 
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_ - 21/2 
2 k a[— cos $ + (l - {~) sin $) } 

f f J p„ ck u e J<wt) 


n j P c ck u 

— TT- 


[l - Jkr] dr d$ 


(A-l) 


The first integral of this is found readily and is 


p Q ck u e^ 0 * 


inf 

K*^ 


U ♦ (—) (cos 2 4 - sin 2 $) 

cL 


2 1/2 2 1/2 
+ 2(—)cos $(1 - (^) sin 2 $) ] + ’a[— cos f + (l • (f-) sin 2 $) ] / d$ 


(A-2) 


The second term in the first bracket integrates to zero because of 
symactry. Hence the pressure is 


o cu e J<at 

Pi = “ 2 — £ t(ka) 2 ♦ E (I 2)] 


(A-3) 


When the elliptic integral is expanded by means of (2k) this becomes 


. e 


[(tat) 2 ♦ 2 } ka {1 - |<2) 2 - J^f) 4 ♦ ...}] 


<A-k> 


This result is exact to the first power of (ka). 

The force on the disk is found by integrating over the surface 


„ /k Jwt 
P 0 cu e° 


cn a 

| | ((ka) 2 + 2J kaU - £(2) 2 _ |j r (|) ,t }]o do d* 


(A-5) 
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where \p is measured about the center of the disk. This can be 
integrated directly to give 


F 


it P Q cu e 
2 







(A-6) 


The imaginary part of this equation can be used to calculate the 
apparent mass. With k = ~ 

F imaginary ' >* “ 8 ^ ^ 


The term in brackets is the acceleration of the disk. Hence the 
apparent mass is 


» a = 2.699 P 0 a 3 (A-8) 

Equation (?) can be transformed to the pendulum form by inserting 
(51a) in place of u. Then equation (A-U) still holds f ~ Pj with 

a 

(i -a sin \|»)6 instead of u and cm additional ten- jcessary to 
c • ^ • 

complete the description. This term is 

„ „ 2 1/2 
2n a[— cos $ + (l -(— ) sin^ $) ] 

* . fit fit 

re‘ Jkr sin (4> - i|»)dr d<|> (A-9) 

0 0 

Expanding the exponential and keeping the first order terms as before 
this equation integrates over r to the form 


p II * - 


J p Q ck 0 e 
inr 


jort 


ii 



6*4 


P II " * 


J p Q ck 0 e 


Jut 2it r 


4TT 


l 


k&- 


((f) 


3 3 a 2 2 

coa $ + 3(—) cos $(1 

ft ft 


2 j^/2 2 

- (f) sin 2 <j>) ♦ 3(f) cos $(1 - (-) sin 2 $) 

ft a a 

2 3/2 2 2 

+ (l - (f) sin 2 $) } + {(—) (cos 2 $ - sin 2 $) 

a c a 

2 

♦ 1 + 2(f) cos $(1 - (-) sin 2 ♦)}} sin (♦ - 4») d+ (A-10) 

a a 

as 'before, making use of the identity- 

sin ($ - 40 = sin $ cos ^ - cos $ sin 4* 


and integrating, most terms drop out; and the result is 

2 


P o c0 e 

p n * 5HT 


jut 


2/P, . , „/* o 


1 * 0 


t(ka)"a(f) ♦ J ftkaja^ 


«f> 


E 'f- r> 


i- (?> 


v — kc|, f )) 


<f> 


sin 4> 


(A-ll) 


With the elliptic integrals expanded by (23) and (2*4) this becomes 


P n c6 e 


Jut 


2 b 

I/O, I/O, 


■II 


T (f)[ (ka) a ♦ j(ka)a(l - f(*) - )) (A-12) 


Pj and p^ are then added to get the total pressure 
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P„c8 


free space 


U c.g. (km > 2 * J<1 “ )U c.g. ,2 -2 < f ) -§$ 


- o gin *{1 - §<f> - |c<f) >1> 


The torque equation is obtained by integrating the pressure 
equation over the disk with the appropriate moment am 


t = 


P Q c0 e 


jut 


2ir a 

J ( " to)£ ‘eg. ♦ J““» l =.g.< 2 - W 


0 0 


2 1 » 

- a sin <j>(l - - |u^“) )]) o sin $ a da d<; 


Integration over <p first gives zero for all terms containing 


l . Thus the result is 

C • £ • 


T = - 


rr p Q c0 e 


jut 


a 

f -J (ka){ ° 3 - 1 H "IfV da 

J a a 


- J w 0 eJwt 


Then dividing by (90) the attached inertia becomes 


I 


a 


0.2792 



4 

> 

(A-13) 


(A-lU) 


(A-15) 


(A-16) 



APPENDIX B 


SOLUTION TO THE EQUATION OF MOTION OF A PENDULUM WITH AERODYNAMIC DRAG 


A perturbation solution to equation (50) is given in this appendix. 
This solution is due to E. M. McDaid and is given he e for the sake of 
completeness. 

Equation (5C) can be written 


9 ♦ 


p(t ) 3 scl 

aero T> 


e|e| 


m g £ 




6*0 


(£-1) 


Nov let 


e 


p(£_ A _ n ) 3s C_ 
aero D 


(B-2) 


and 


C 


2 


m 

s 



(B-3) 


Then equation (B-l) becomes 

e + e e | e | + c 2 e * o o < t « l (b-4) 

where the fact that 0 < e « 1 can be verified by the substitution 
of numerical values. 

Extend the domain of the function to the additional variable 
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t » e v g(t) 


(B-5) 


How apply the principle of maximal balance to determine v: 


d0 . li . -V ' ae 
dt ~ 3t e 8 3 t 


(B-6) 


d 2 0 _ 3 2 0 


v • a 2 © 


a “" 2c g ltl ; r + e "' «" tt + e ' 6 ST 


2 2 
dt* at* 


_2v *2 3^6 . _v - 30 
2 


(B-7) 


3t 


Then applying these derivatives back in the original differential 
equation results in 


3 2 a v * a 2 e 2 v * a 2 e 

' + 2e g + € g ■ ■ ~ + e 
3t 2 6 3t3t R S T 2 


v l 30 . c /30\ 2 . pr v+l * 30 30 
g aPT + 6 \3t/ + 2e 8 3t 3? 


. _2v 

+ e g 


*2 /30\ 
8 r) 


+ 6 2 0 


°-H >0 

(B-8) 


or 


is + 2t « ; is. + c 2v • is. + e v ~ a 

at 2 * E g 3t3x E g 3t 2 e 8 3t 


_2\H-1 

- e g 


- e l%\ 

w 

*2 /30\ 

8 W 


. 2e ^i ;ii m 

* 6 8 3t 3 t 


+ 6 2 0 = 0, || < 0 


(B-9) 


Consider the algebraic equations. They contain terms in the 
following powers of £: e* 3 , e v , e 2 ^, e, e v+3 ‘, and e 2 ^^. Solving 


for v by equating exponents of varioiis pairs of terms, the possible 



values of v are: -1, - j, 0, and 1. The c terms thus 

become, for each choice of v: 


B 

1 

It 


V = 1 


0 

0 

0 

0 

c o 

e 

e 

£ 

e 

£ 

c° 

-i 

e 

i 

e " 2 

£ 

i 

e" 2 

0 

•2 

l 

2 


e 

e 

e 

£ 

e 

1 

e 

£ 

e 

£ 

1 

e 

1 

0 

1 

2 

l 

e 

e 

e 2 

e 

e 

e 1 

e- 1 

£ 

e 3 

1 

e 2 


The best balance is yielded by v = 1 (the case v * 0) is the trivial 
case of no scaling at all). Thus (B-5) becomes 


x * e g(t) 


Nov assume a solution of the form 


0 * 0 t (t, tje 1 

J-0 


Then 


(B-10) 


(B-ll) 
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St * It * e « If * 9 ot * £(e n * « V * e2(9 a * « V * ••• <B - 12) 


If 9 » °» c - e< 9 ot > 2 - e 2 < 29 ~< 9 i. + « 9 „>> * ••• 


'Ot' v lt • B W 1T' 


(B-13) 


If 9 ‘ E H# ■ - £ ( 9 0t) 2 - * 6 ... <B-)M 


at' It ° Or' 


0 * 0 * 29 9 & * £ 9 H *■* * 2 J * 9 ott * £ < 9 itt ♦ *<W ♦ 9 v 


+ e ^?tt + 2g6 ltT + *®lT + g e OXT^ 


(B-15) 


Thus if 0 > 0 


(6 Ott + * e[6 ltt + + 2g6 0tt + « 0 OT * *lt ] * e2( *** ) + * 3( — ) 


+ . • * * 0 (B-l6) 


The resultant system to second order is 


6 0tt + 5 6 0 * ° 


(B-17) 


W “ t2g8 0tT + »<* * W 


(B-l8) 


with the initial conditions at t * 0 
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V° • 0 ot * e o 


(B-1S) 


* 0 


"it * * 6 0T * 0 


Also, if 0 < 0 


(9 0tt * ^o> * £(e itt * * 2e9 0tx * « e ot 


e2 0t> * 


) + 


+ ... =0 (B-20) 


This system of equations is 


6 0tt + ‘ 6 0 = ° 


e itt + - - [2g0 otT ♦ - e^] 


(B-21) 

(B-22) 


with the initial conditions at t«0 


e 0 ■ 0 • e ot = 6 o 


(B-23) 


e x * o 


0 + c0 * 0 

it * on 


The above two systems of equations are equivalent to 
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(B-24) 


Vt * “ 0 




9 xtt * ' 2 «1 - 

^ 2g0 Otx + g e ox * l 0 oJ 0 ot^ 

(B-25) 

with the initial conditions 





O 

1! 

O 

CD 

e ot * *o 

(B-26) 



6 i s ° 

* e it * ge ot * 0 


Solving for 

e o 






e Q » A q (x) sin(ct) 

(B-27) 

The function 

A o 

is determined 

from the second equation: 



6 ltt * 0 1 * ~ [2«^<T)co8(ct) + gA^T)sin(ct) + A^T^cosUt ) |cos(r,t ) | ] 

(B-20) 


In order for secular terms to be absent, 

2tt 

sin(Ct){2gCA^(x)cos(ct) + gA* (x)sin(ct) 

0 

+ Aq( x ) £ 2 cob ( £ t ) | cos ( £t ) | }dt = 0 (B-29) 


T1 



and 


2n 

j" cos ( £t ) { 2g£A^ ( x ) cos ( *t ) + gA'(x)sin(5t) 

0 

+ A 2 (x)c 2 cos(ct) |cos(£t) | }dt = 0 


First, choose g - 0 

This implies g = £ since the integration constant is arbitrary 


8 ■ Ct 


and 


r = e?t 

This simplifies the above equations to 

2n 

j sin(;t){2C 2 A^(T)cos(ct) + S 2 A 2 (x)cosUt) |cos(ct ) | }dt * 0 
0 

27T 

j co8(tt){2C 2 A^(x)cos(ct) + c 2 AQ(x)eos(ct)|cos(ct)|}dt * 0 
0 

The first equation is satisfied regardless of the value of A^. 
second equation yields 


(B-30) 

(B-32) 

(B-33) 

(B-3U) 

(B— 35 ) 
The 


72 



(B-36) 


2itc 2 A^(t) + | C 2 A 2 (t) = 0 


The frequency term drops out. Hence 


A o< T > 

A=(t) 


u_ 

3* 


Thus 




~ ^ K 
3 ^ 


Uect 


3tr 


+ K 


or 


H 


0 


t) 



It is now necessary to solve for K. From (B-27), taking the tir 
derivative 


38- 3A-(t) 

* 6A 0 (x)cos(Ct) + — ^ — sin(u) 

If this is evaluated at t “ 0 


(3-37) 


(B-38) 


(B-39) 


(B-1*0) 
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**! t»o * tA ° <0> 


But from the initial conditions (B-26) this is equal to 0 Q , so 


V 0) = f 


also from (B-j9) at t = 0 


V c > ■ f 


Hence 


K-f- 

0- 


So, in (B— 39) 


V T> * r — 7 k 


m 


The final solution is then 


6 .Si 

0 c 


A- y ain(Ct) * 0(e) 

1 * HH 


It ran be seen that the amplitude is strongly influenced by the 


frequency. However, the period remains essentially th* same as for 


7 ^ 



undamped motion. This solution is essentially the same as that of 
Bogoliubov and Mitropolsky (ref. 15) , p. 75-77- However, it is thought 
that the method used here is somewhat simpler and more elegant. 
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